Theory of inelastic scattering from magnetic impurities 
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We use the numerical renormalization group method to calculate the single particle matrix ele- 
ments T of the many body T-matrix of the conduction electrons scattered by a magnetic impurity 
at T — temperature. Since T determines both the total and the elastic, spin-diagonal scattering 
cross sections, we are able to compute the full energy-, spin- and magnetic field dependence of the 
inelastic scattering cross section, <7i no i(tL>). We find an almost linear frequency dependence of o"i ne i(w) 



below the Kondo temperature, Tk, which crosses over to a 



behavior only at extremely low 



energies. Our method can be generalized to other quantum impurity models. 
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Quantum mechanical phase coherence in mesoscopic 
structures is destroyed due to inelastic processes, where 
excitations such as spin waves, electron-hole excitations, 
phonons, etc., are created in the environment with a cer- 
tain probability, thus leading to dephasing and loss of 
quantum coherence after a time ~ t v pj. In some weak 
localization measurements of the dephasing time t v down 
to very low temperatures, a surprising saturation of t v 
has been observed Q. This unexpected saturation re- 
mained a puzzle for a long time until recently, when 
further experiments on mesoscopic quantum wires con- 
firmed that the most likely candidates to produce this 
surprising behavior are magnetic impurities 0, 0] . These 
magnetic impurities seem to be present even in samples 
of extreme purity, and unavoidably lead to inelastic scat- 
tering and the dephasing of charge carriers. 

Theoretical calculations also confirmed these expec- 
tations and showed that the experimental data can be 
quantitatively explained assuming weak inelastic scat- 
tering off Kondo impurities 0, 0|. These calculations, 
though, were performed in the weak coupling regime, 
i.e., at energies higher than the Kondo temperature, Tk- 
However, Nozieres Fermi liquid theory teaches us that 
well below Tk the magnetic impurity spin is screened 
by the conduction electrons, and there it acts simply as 
a strong but conventional potential scatterer, and thus 
produces no inelastic scattering. Therefore the inelastic 
scattering rate from magnetic impurities is expected to 
show a peak around Tk and then drop to zero well below 
T K U 

These observations motivate us to study the complete 
energy dependence of the inelastic scattering rate off a 
magnetic impurity. Here we shall focus on the simplest 
possible case of T = temperature, where the inelastic 
scattering rate can be defined as follows: Assume that we 
have a single scattering impurity at the origin and we cre- 
ate an incoming flux of electrons with momentum k, spin 
a, and energy E above the Fermi energy, far away from 
the origin. This incoming flux can be scattered off the 
impurity in two different ways: (i) Either the electrons 
scatter elastically (both energy and spin unchanged) with 



a scattering cross section <j e \(E) into an outgoing sin- 
gle particle state, without perturbing the environment, 
or (ii) they scatter off inelastically with a corresponding 
cross section er; nc i(i?), i.e., and they leave behind some 
electron-hole or spin excitations. 

In the present paper, we show how the full energy and 
magnetic field dependence of <7i no i (E) can be determined. 
The basic idea is simple: The single particle matrix ele- 
ments of the many-body T-matrix, (k <r|T|k' a'), deter- 
mine the elastic cross section, but they are also related 
to the total scattering cross section, a tot — c e i + <7i ne i 
through the optical theorem. Therefore, we only have to 
find a way to compute the (k c|T|k' er'}'s to obtain the 
inelastic scattering cross section as the difference of the 
total and elastic scattering cross sections: 

°tael = °total ~ a c\ ■ (!) 

To determine (k <j|T|k' cr'), we shall first relate them 
through the reduction formulas to some local correla- 
tion functions ||, which we shall then calculate using 
the non-perturbative method of the numerical renormal- 
ization group (NRG) 0. Though here we shall focus 
exclusively on the case of T = temperature, where ex- 
citations are created from the vacuum state [lfj , our dis- 
cussions carry over, with some modifications, to the case 
of finite temperatures too [llj| . 

To be specific, we first consider the Anderson model, 
but our method is rather general and applies to practi- 
cally any local quantum impurity problem. We write the 
Hamiltonian as H = Hq + H lnt , where Ho denotes the 
'free' quadratic part of the Hamiltonian, 

Ha = Y^ ZddUa + ^ f (k)4 a c k a , 
a <T,k 

and H m t stands for the on-site Hubbard interaction and 
hybridization 

H int = U n T 7H + VJ2( C 1 <A + h - c -) » ( 2 ) 



The operator creates an electron 



with n„ = d^do 
in a plane wave state with momentum k, energy £(k) = 



£- — /i and spin a, while d a is the annihilation operator 
of the d-electron. 

We proceed to define incoming and outgoing scatter- 
ing states as well as the corresponding field operators and 
Hilbert spaces 8]. As the impurity is local, the interac- 
tion switches off far away and the 'in' and 'out' states 
are eigenstates of the full Hamiltonian with the asymp- 
totic boundary condition of behaving as plane waves in 
the x — y — oo and x — > oo limits, respectively. The many- 
body S-matrix and the T-matrix elements are then sim- 
ply defined in terms of the overlaps of the incoming and 
outgoing scattering states, 



(fr, out | a, in) 
§=l+if 



(b, w\S\a, in) 



(3) 
(4) 



In the interaction representation, the explicit form of 
the S-matrix is given by the well-known expression S = 
Texp{— i J_ iJjnt(i) dt}, where T is the time ordering 
operator. 

Since we are primarily interested in the single-particle 
matrix elements of T, we consider the case where asymp- 
totically | a, in) is simply a single particle low energy exci- 
tation of the free vacuum defined by H (x — — oo) — > Hq, 
having momentum k and spin a: |o, in) = |k, a) . Then 
the single-particle matrix elements of the T-matrix read 

<k, *\f\k>, a') = 2^(C(k) - C(k'))(k, a\T\k\ a') , (5) 

where we separated a Dirac delta contribution due to 
energy conservation and defined the on-shell T-matrix 
T. Next, following the manipulations of Ref. we re- 
express the off-diagonal (k ^ k') matrix elements of the 
T-matrix as: 



(ka|T|kV) = 



(0) 



<(£,sk )S k') GoH^sk') 



where s = sgn £ distinguishes electron-like excitations 
from hole- like excitations, Gq 1 = ijj + fi + de- 
notes the inverse of the free Green's function, and G is 
the time-ordered single particle Green's function. The 
meaning of Eq. (j^J becomes more transparent in the dia- 
grammatic language of Fig.^ As indicated by the large, 
thin crosses there, one has to drop the contributions of 
the two external legs of all scattering diagrams to the 
single electron Green's function and the rest is just the 
on-shell single particle matrix element of the many-body 
T-matrix. In the particular case of the Anderson model, 
T does not depend on the direction of incoming and out- 
going momenta, and a simple Dyson equation relates it 
to the d-level's time-ordered propagator (see Fig. ^) 



T a (cu) = -sV 2 



Gs-cr 
d 



where s = sgn u> and we allowed for spin-dependence due 
to the presence of an external magnetic field B. Eq. Q 
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FIG. 1: (a) Diagrammatic derivation of Eq. Dashed and 
continuous lines denote the bare propagators of the d- level 
and the conduction electrons, small fat crosses stand for hy- 
bridization V, and wavy lines denote the on-site interaction 
U. (b) Diagrammatic representation of the T-matrix in the 
Kondo problem. Dashed lines denote pseudofermion propa- 
gators and describe the evolution of the impurity spin, while 
continuous lines denote free conduction electron propagators. 
Filled circles stand for the exchange interaction J. The first 
term of the T-matrix is simply proportional to the expecta- 
tion value of the impurity spin, it is frequency independent, 
and vanishes in the absence of magnetic field. The second 
term can be identified as the composite fermions correlation 
function. 



has also been derived in a somewhat different way in 
Ref. [H 

According to the optical theorem, the spin-dependent 
total scattering cross section is given by the imaginary 
part of the diagonal matrix elements of the T-matrix: 



"toted = — Im(pcr|T|pcr) 

V F 



(8) 



where vf denotes the Fermi velocity. The elastic scat- 
tering cross section, on the other hand, is related to the 
square of T: 



05 = $y 3 ^W-0\(p'*\T\p<7)f 



(9) 



Once these two cross-sections are known, we can compute 
the inelastic cross section <Ji ne \ through Eq. QJ. 
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It is instructive to rewrite cxinei in case of the Anderson 
model. For electrons we have 



4tt 



incl 



(u > 0) = vo- 



k' 2 

F 



IGSI 2 



(10) 



where T = 2irV 2 go denotes the width of the d-level and 
Qa = kp/2TT 2 VF is the conduction electrons' density of 
states for one spin direction. For B 
reduces to 



0, this expression 



ei(w > 0) 



2tt 



r (-£») 



k 2 F (u 



e d -E'(u;)) 2 + (£"M-§) 2 



where £' and S" denote the real and imaginary parts of 
the d-propagator's self-energy. The analytical properties 
of the Green's function imply that the above expression 
is always positive and it only vanishes where S" becomes 
zero. Furthermore, the Fermi liquid theory of Yamada 
and Yoshida tells us that S" ~ to 1 as to — > 0, and thus 
the inelastic scattering rate vanishes as io 2 at the Fermi 
energy. Note that at the same time the total scattering 
cross section approaches the unitary limit. 

To compute the full behavior of cri ne i we determined 
the T-matrix using the numerical renormalization group 
method Since we were dominantly interested in the 
low-energy universal regime of the Anderson model, we 
took the U/2 = — — > oo limit and performed the cal- 
culations using the Kondo Hamiltonian, 



(11) 



k.k' 



Up to an overall normalization factor, the low-energy 
part of the spectral function of the original d-level prop- 
agator in the Anderson model can be shown to corre- 
spond to the spectral function of the following compos- 
ite fermion operator in the Kondo model [l2(, F a = 
J2 a ' k S-daa'Cka' . (For a diagrammatic proof, see Fig.^l. 
The imaginary part of G ( i ~ T can thus be deter- 
mined by simply computing the spectral function Qf(u) 
of the composite fermion numerically, and then a Hilbert 
transform can be used to get the real part of Gd and 
thus the full T-matrix. In all these calculations it is 
essential to have high quality data. It is also crucial 
to determine the normalization factor of TGd correctly. 
This can be done by using the Fermi liquid relation, 



-Im TG a d (uo = 0+) = 2shr S a , with S a the phase shift at 
the Fermi energy. We extracted the latter directly from 
the finite size NRG spectrum of the Kondo model [HEl- 
Our results for the case of no external magnetic field 
are shown in Fig. [21 Most of the scattering is inelas- 
tic at energies above the Kondo energy, \w\ > Tk- De- 
creasing the energy of the incoming electrons (holes), 
the total scattering cross section increases and, at en- 
ergies below the Kondo scale, it finally saturates at a 
value do = 4:ir/k F . This behavior must be contrasted 
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FIG. 2: Inelastic, elastic, and total scattering cross sections 
in units of ao = Air /k^ at T — and in the absence of mag- 
netic field, as a function of the logarithm of the energy of the 
incoming electron. Only the electronic contribution (ui > 0) 
is plotted. (Ti nc i only has a very weak (logarithmic) energy 
dependence above Tk, scales approximately linearly with uo 
for 0.01 T K <uj <T K , and scales as ~ lu 2 for lu < 0.01 Tk- 
The inset shows the ~ u> 2 and <Ti nc i ~ 1/ln 2 {Tk/u) regimes 
for iv Tk and u> 3> Tk, respectively. 



to the inelastic scattering rate, which slowly increases 
as u) decreases, has a broad maximum around Tk then 
suddenly drops and vanishes at the Fermi energy. On 
linear energy scales (see Fig.|3J), Omd varies rather slowly 
above Tk, is very large even at to ~ 2QTk, and vanishes 
rather suddenly around u) ~ Tk- For very small energies 
tTinci ~ uj 2 , in agreement with Fermi liquid theory, how- 
ever, this quadratic behavior appears only at very low 
energies, and <7i no i is almost linear for 0.1 Tk < oj < Tk- 
At energies to 3> Tk the inelastic rate is simply dom- 
inated by energy-conserving spin-flip scattering, and is 
therefore expected to scale as ~ l/ln 2 (Tft-/cc>), as we in- 
deed find numerically. Note that the Nagaoka-Suhl ap- 
proximation 0| is only appropriate for u> S> Tk (see 
inset of Fig. 2). 

We also computed the inelastic scattering rate in the 
presence of a local magnetic field, B, directed downwards 
along the z-axis (see Fig.|3J). In this case there is a dra- 
matic difference between the inelastic scattering proper- 
ties of spin up and spin down particles. Already a small 
field, B ~ O.ITk results in a strong spin-dependence of 
the inelastic scattering, but for B ~ Tk, this difference 
is even more dramatic. At this field the spin of the mag- 
netic impurity is practically aligned with the external 
field and points downwards. Therefore an incoming spin 
down particle (electron or hole) is unable to flip the im- 
purity spin. More precisely, only higher order inelastic 
processes can result in a flip of the local impurity spin. 
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FIG. 3: Energy dependence of spin-dependent elastic and 
inelastic scattering rates in units of o~q = itr/kp, at T = 
and in the presence of a local magnetic field B. 



This is, however, not true for spin up particles: An in- 
coming spin up electron can exchange its spin with the 
magnetic impurity in a first order process, resulting in a 
maximum in the inelastic scattering cross section around 
u> ~ B for spin up electrons and holes and a very broad 
inelastic background for w > B. 

Though here we mostly focused on the simplest cases 
of the Anderson and the single channel spin S = 1/2 
Kondo models at T = temperature, our formalism can 
lie extended to other models and to finite temperatures as 
well [l"H ]. I n particular, while for some quantum impurity 
models no simple diagrammatic theory is available, the 
composite fermion's spectral function can be computed 
in any Kondo-type model to obtain the matrix elements 
of the T-matrix, and the renormalization group flow of 
the ei genv alues of the S-matrix can be studied in all these 
cases jllj. While usually a thorough numerical analysis 
is needed to understand the full behavior of <7i ne i, in some 
models simple analytical results can also be obtained. In 
the specific case of the 2-channel Kondo problem, e.g., 
we know that the single particle matrix elements of the 
S'-matrix identically vanish at the Fermi energy, lu = 
[3 El This implies that g T 2CK {u = 0) = i/n and 
leads to the rather surprising relation at the Fermi level, 



T 2CK 



T 2CK 



°fot /2 ; Though the S'-matrix vanishes 
identically, half of the scattering processes remain elastic. 



The non-vanishing inelastic scattering rate is character- 
istic of non- Fermi liquid quantum impurity models. Ap- 
plication of any finite magnetic field drives the 2-channcl 
Kondo model to a Fermi liquid fixed point, and gives 
rise to a vanishing inelastic scattering rate at the Fermi 
energy. 

We have to emphasize that, though they must be re- 
lated, the inelastic scattering rate we computed is not 
identical to the dephasing rate measured in weak lo- 
calization experiments 0, since the former incorporates 
spin-flip scattering processes as well as the creation of 
electron-hole pairs. While we only computed <7i nc i(w, T — 
0), we expect that <7i n ei(^ = 0, T) has a very similar form. 
In this sense, our finding that the inelastic scattering 
rate is roughly linear with uj for 0.05 Tk < oj < 0.5Tr- 
agrees qualitatively with the recent experimental results 
of Ref. [H 

In summary, we have shown how the full energy and 
magnetic field behavior of the T — inelastic scatter- 
ing rate can be computed by exploiting the reduction 
formulas and then using the powerful machinery of nu- 
merical renormalization group to compute the single par- 
ticle matrix elements of the many-body T-matrix. We 
have shown that the Fermi liquid theory of Yamada and 
Yoshida directly implies a quadratically vanishing inelas- 
tic scattering rate at the Fermi energy in the specific 
case of the Anderson model. Scattering properties of 
the Kondo model have been computed by calculating the 
composite Fermion's spectral function. Our numerical 
calculations show that the abovemcntioned cr inc i ~ uj 2 
regime appears only at energies well below the Kondo 
temperature. In a magnetic field B > Tk the inelastic 
scattering is very sensitive to the spin of the scattered 
single-particle excitation. 
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